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ABSTRACT 

Context. Intermediate and high velocity HI clouds rain onto the plane of our Galaxy. They are observed at heights of 
between 500 and 1500 pc, falling onto the Galactic plane at velocities from 50 to 140 km s~^. 

Aims. To explain the origin of these clouds, we present a galactic fountain model, driven by the wind from a super 
stellar cluster (SSC). 

Methods. We solve the equations for a steady, radiative de Laval nozzle flow. We consider two effects not considered 
previously in astrophysical nozzle flow models: cooling functions for different metallicities, and the direct action of the 
galactic gravitational field on the gas flowing along the nozzle. 

Results. For an adiabatic nozzle flow, the gravity acting directly on the gas within the nozzle "stalls" the nozzle flow for 
initial wind velocities lower than the escape velocity from the Galaxy. For the same wind velocity, a radiative nozzle 
flow stalls at lower altitudes above the galactic plane. We find that SSC winds with velocities of — 500 — 800 km s~^ 
produce nozzles stall at heights of Xm = 1 — 15 kpc. The stalled nozzle flow then rains back onto the galactic plane at 
velocities in the range observed in intermediate and high velocity HI clouds. 

Conclusions. We study a nozzle flow driven by a wind from a SSC close to the Galactic centre. We flnd that for velocities 
within the range expected for a SSC wind, we can produce nozzle flows that stall above the galactic plane. These stalled 
flows produce cool, infalling clouds with velocities similar to those of intermediate and high velocity HI clouds. 

Key words. Hydrodynamics - ISM: jets and outflows - Galaxies: star clusters - Galaxies: starburst - Galaxies: intergalactic 
medium 



1. Introduction 

The observations of Munch & Zirin (1961) showed absorp- 
tion (mostly Ca II) lines of interstellar clouds at distances of 
500-1500 pc above the Galactic plane. These high altitude 
interstellar clouds have been described by many authors 
as neutral hydrogen (HI) clouds at distances of between 
300 pc and 12 kpc with infall velocities between —50 and 
-140 km s-\ 

The clouds have been divided into two categories: 
intermediate-velocity clouds (IVCs) and high-velocity 
clouds (HVCs). IVCs have radial velocities in the range of 
30 km s~^ < |visi^| < 90 km s^^ with respect to the LSR, 
and HVCs have radial velocities |vis/j| > 90 km s~^. In the 
Galaxy (and in other galaxies), the IVCs have typical alti- 
tudes of between 300 pc and 2.5 kpc, and masses of up to 
~ 10^ Mq. HVCs have been observed at higher altitudes (of 
up to ^ 12 kpc). Several complexes of HVCs with masses 
of ~ lO'^ — 10^ M0 have been observed (see, e. g., Wakker & 
van Woerden 1997) . The general characteristics of HI clouds 
in the Galaxy were discussed by Dickey & Lockman (1990) 
and Wakker & van Woerden (1997). Their general chemical 
properties were described by Lu et al. (1998), Wakker et al. 
(1999), Murphy et al. (2000), Sembach et al. (2000), and 
Bluhm et al. (2001). 

Significant progress has been made in exploring the dis- 
tribution of IVCs and HVCs in the halo of the Galaxy (see, 
e.g., Dettmar 2005). However, the origin of these clouds is 
unclear. Oort (1970) proposed that high altitude clouds (in 
the Galaxy) are formed in past episodes of star formation. 



Several authors have used this scenario to explain the origin 
of the IVCs. However, the HVCs have been considered to 
be extragalactic objects that are about to merge with the 
Galaxy (Fraternah & Binney 2006, Mehoh et al. 2008). An 
alternative scenario was proposed by Blitz et al. (1999), 
who suggested that HVC and IVC could be relics of the 
early stages of Galactic formation. 

It has also been proposed that both IVCs and HVCs are 
formed in the disk of our galaxy, in the so-called "galac- 
tic fountains" . These models were described by Shapiro & 
Field (1976), Houck & Bregman (1990), and de Gouveia 
Dal Pino et al. (2009). In this scenario, the cloud formation 
is fed by chimneys of hot, ascending gas that cools and falls 
back onto the Galactic disk in the form of discrete clouds. 
Observations indicate that chimneys are indeed generated 
in other galaxies by multiple supernovae (SNe) explosions 
(Veilleux, Cecil & Bland-Hawthorn 2005, Konstantopoulos 
et al. 2008). 

De Gouveia Dal Pino et al. (2009) showed that SNe 
explosions can lead to the formation of hot superbubbles 
that drive a supersonic wind out of the galactic disk. From 
3D, radiative hydrodynamical simulations, they found that 
the galactic fountains driven by SNe may reach altitudes 
of up to ~ 5 kpc in the halo and thus allow the formation 
of IVCs. However, the HVCs (found at higher altitudes) 
require a different production mechanism. 

We propose that HVCs (and possibly also some of the 
IVCs) could be formed by fountains driven by super stellar 
clusters (SSCs). SSCs are dense clusters of young massive 
stars. These stellar clusters contain hundreds or thousands 
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of very young, energetic stars, and have stellar densities far 
greater than those seen in normal OB associations. These 
star clusters have ages of ~ 1 — 10 Myr, radii of ^ 1 - 
10 pc, total cluster masses of 10^ - 10^ Mq (Melo et al. 
2005 reported a mean mass per stellar cluster of 2 x 
10^ Mq, for M82), and typical cluster wind velocities of 
500 - 1500 km s'^ range (see also, Silich et al. 2004). The 
central stellar densities of SSCs have values of up to ~ 10^ 
M0 pc~^. However, SSCs with older ages and/or higher 
masses do exist (Walcher et al. 2006 reports a cluster with 
6 X 10^ Mq). On the other hand, the metallicity of the gas 
reinserted (via supernova explosions and stellar winds) into 
the stellar cluster medium, reaches supersolar values (1 < 
Z < 14, see also Meynet & Maeder 2002 and Tenorio-Tagle, 
Silich, and Mufioz-Tufion 2003) during most of the cluster 
evolution. 

SSCs have been observed in a wide range of star-forming 
galaxies, such as merging systems (NGC 4038/4039, 
Whitmore & Schweizer 1995), dwarf galaxies (Henize 2-10, 
Johnson et al. 2000), classical starbursts (M82, Gallagher 
& Smith 1999, Melo et al. 2005), as well as in the Galaxy 
(Arches Cluster: Nagata et al. 1995; Cotera et al. 1996; 
Serabyn, Shupe & Figer 1998) amongst many other sys- 
tems (for a review, see Whitmore 2000). 

In this paper, we present a model for the formation of 
low and high velocity HI clouds by means of radiative noz- 
zle flows driven by the wind from a SSC. The nozzle flow is 
produced by the interaction of the shocked SSC wind with 
the stratified halo and the gravitational potential of the 
Galaxy. In our model, we assume that the SSC that pro- 
duces to the fountain is close to the Galactic centre. This 
is an appropriate assumption if we consider the wind from 
the Arches cluster or from other possible massive (but still 
undetected) stellar clusters close to the Galactic centre. In 
galaxies such as M82 and NGC253, it appears that we do 
see outflows ejected from SSCs in the central regions of 
these galaxies (see Rodriguez-Gonzalez et al. 2008). 

The paper is organized as follows. In Section 2, we 
describe the equations for a radiative nozzle flow, and in 
Section 3 we discuss the analytic solution obtained in the 
adiabatic case. In Section 4, we describe the Galactic po- 
tential and hot gas distribution, which are used in Section 
5 to compute numerical integrations of the radiative nozzle 
flow equations. Finally, in Section 6 we discuss the appli- 
cation to these models for low and high velocity HI cloud 
formation, and we present our conclusions in Section 7. 

2. Equations for a radiative nozzle flow 

We assume that the interaction of a wind from a SSC (lo- 
cated close to the centre of the galaxy) with the surround- 
ing interstellar medium (ISM) leads to the formation of 
an approximately spherical shock. The thermalized, post- 
shock material is then channeled into two oppositely di- 
rected nozzles because of the pressure stratification of the 
environmental material. In the context of the generation of 
jets by active galactic nuclei, this scenario was first explored 
by Blandford & Rees (1974) 

We consider a one-dimensional, stationary, radiative 
nozzle flow under the action of a gravitational potential. 
We first write the energy equation in the "entropy conser- 
vation" form 



dx 



Pp- 



(7 - l)L 



(1) 



where x is the direction of propagation of the flow, P is 
the gas pressure, p the density, u the flow velocity, 7 the 
specific heat ratio, and L is the energy loss per unit time 
and volume. We now consider the momentum equation in 
the form 



1 dP 
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where ^{x) is the gravitational potential, which we take to 
be a known function. Equations ([1]) and ^ can be com- 
bined to obtain 
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Equation ([3]) simplifies to Bernoulli's theorem for the case 
of a non-radiative fiow (i. e., for L = 0). 

Following Blandford & Rees (1974), we now assume that 
the nozzle fiow has a slowly varying cross-section, so that 
the pressure of the flow always adjusts to a value close to the 
pressure Pe{x) of the surrounding environment. If one also 
knows the cooling function L as a function of the density 
and temperature of the flow, the nozzle flow can then be 
obtained by integrating the system described by eqs. ([1]) 
and (O. 

The initial conditions for this integration are obtained 
by assuming that the initially isotropic wind from the clus- 
ter (of mass-loss rate M and terminal velocity Vw) follows 
a strong, spherical shock before flowing along the nozzle. 
Assuming that the shock is strong, we have a post-shock 
velocity 



7+1 



and a post-shock sound speed 

_ v/27(7 - 1) 
7+1 



(5) 



(6) 



Following Raga & Canto (1989), we assume that the initial 
velocity and sound speed along the nozzle flow are vq « u„.. 



and Co 



respectively (see Eqs. [5] and [6|). The initial 



density po — p{x — 0) of the nozzle flow can be obtained 
from the P{x) = Pe{x) lateral pressure-balance condition 
(see above) 



Po 



(7) 



ps 



where ^^6(0) is the (known) environmental pressure ai x = 
0, and Cps is given by Eq. ([6]). The initial conditions for the 
integration of the nozzle flow problem (Eqs. [1] and [3]) are 
then given by Eqs. (OE]). 

Finally, the mass conservation equation for a flow of 
slowly varying cross-section 



M = TTRyu , 



(8) 



where M is the mass-loss rate, can be used to calculate the 
cylindrical radius R of the cross-section as a function of 
position X along the nozzle flow. 
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Fig. 1. Free-fall velocity as a function of Xm for different 
final heights xq — 0.3, 2.5, 5, 10 and 15 kpc (solid, dotted, 
dashed, dash-dot, and dash-dot-dot lines, respectively). 

3. The adiabatic solution 

For a non-radiative flow (i.e., with L = 0), Eq. ([3|) can be 
trivially integrated to obtain 



u 

T 



7 ■ 



(9) 



where $o = ~ 0), and the constant on the right is 
obtained by setting u = Vpg and = Cps (see Eqs. [5] and 
m) at X 0. 

For L = 0, from Eqs. ([9]), ([2]), dH) and dH) and the lateral 
pressure balance condition, we obtain the solution: 



u\x) 



47 



(7 + 1)^ 



Peix) 



Pq 



(7-l)/7~ 



-2 [$o - Hx)] 



where Pe{x) is the environmental pressure stratification, 
and $(a;) is the gravitational potential. If these two func- 
tions are known, Eq. pU|) directly gives the velocity along 
the nozzle flow as a function of position x. We should note 
that this solution is equivalent to the one of Blandford & 
Rees (1974), but also includes the direct effect of the grav- 
itational force on the material within the nozzle flow. 

It is easy to check that for a; = 0, Eq. pI7|) correctly gives 



uiO) 



(see Eq. [S]), the imposed initial value for the 



flow velocity. For cc — > 00, for all finite mass environmental 
stratifications we have Pe{x) and ^(x) (setting 
the zero of the potential at infinity) . We then obtain 



u{x) — > Moo 

Because $0 
tion only if 



\/vl, + 2$o , ioY x ^ 00 . 



(11) 



< 0, this equation clearly has a physical solu- 



(12) 



cluster of massive stars (with ^ 1000 km s ^) satisfies 



the condition shown in Eq. (jl2p , and therefore produces an 
adiabatic nozzle flow that escapes from the Galaxy. In the 
following sections, we show that this is not necessarily the 
case for radiative nozzle flows. 



4. The Galactic potential and hot gas distribution 

We assume that the halo has a uniform temperature. An 
isothermal halo is an assumption that is consistent with the 
observed diffuse, soft X-ray background emision (Garmire 
et al. 1992). The Galactic gas density pg and thermal pres- 
sure Pg are given by (Espresate et al. 2002) 



30 p,(x)=p,,oexp[*o-*(-)l/s 



and, 

Pg(x) =Foexp[*°-*(")l/=' 



(13) 



(14) 



where, Pq — PgflCg and pg^ are the pressure and density, 
(respectively), of the hot gas at the Galactic centre, and Cg 
is the (uniform) isothermal sound speed of the halo. For the 
gravitational potential ^{x), we use the model proposed 
by Allen & Santillan (1991) and Espresate et al. (2002). 
These authors defined three components for the Galactic 
mass: (1) a nucleus, represented by a spherically symmetric 
distribution, (2) a disk component, (3) a spherical, massive 
halo. Hence, the gravitational potential is given by the sum 
of these three contributions. 

The units are such that if x is in kiloparsecs and Mi, 
M2, and M3 are in Galactic mass units (gmu=2.32x 10^ 
Mq), the units for ^{x) are 100 km^ s~^ (Espresate 2002). 
The equations for the potentials along the rotation axis of 
the Galaxy, for the different components are as follows 
1. the nucleus is represented by 
(10) 

Ml 



where, Mi=606 gmu and bi = 0.3873 kpc; 
2. the disk component is represented by 



$2(X) = 



Mo 



as + (x2-f 62)1/2' 



(16) 



where, M2 — 3690.0 gmu, 02 — 5.3178 kpc, and 62 = 0.25 
kpc; and, 

3. finally, a halo truncated at 100 kpc (see also Allen & 
Santillan 1991), has a potential 



$3(2;) 



1 + ln 



C3, 



(17) 



where, 03 — 12.0 kpc, 63 = 377.04 gmu/kpc, and C3 = 1201 
gmu/kpc. The total gravitational potential at any point x 
along the rotation axis of the galaxy is then given by 



in other words, the nozzle flow extends out to infinity only 
if the velocity of the wind (which is producing the nozzle 
flow) exceeds the escape velocity. If this condition is not 
met, the nozzle flow reaches a maximum height Xm (which 
can be computed from Eq. [TO] setting u = 0), at which it 
stalls. However, since the escape velocity from the centre of 
our galaxy has a value Vesc ~ 700 km s"""^, a wind from a 



$(x) = $i(x) -I- $2(2;) + ^•3(a;)- 



(18) 



This expression for the gravitational potential (from 
Espresate et al. 2002) implies a total mass for the Galaxy 
of Mt = 9.0 X 10"^^ M0, and a local circular velocity of 220 
km s^^ and tq = 8.5 kpc. With this potential, we can com- 
pute the free-fall velocity of material injected at a height 
Xm above the galactic plane (by a nozzle flow originating 
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close to the galactic centre). By the time that the material 
injected (with a velocity « 0) at Xm has fallen to a height 
xo above the galactic plane, it has acquired a velocity: 

vj = ^xm) - Hxo) . (19) 

Figure 1 shows as a function of Xm for different final 
heights Xo, for values of Xo that are typical of the observed 
altitudes of IVCs and HVCs. The resulting values of are 
used in the following sections to determine the conditions 
that have to be satisfied by nozzle flows to produce clouds 
with velocities similar to those of IVCs and HVCs. 

5. The models 

We consider a nozzle flow fed by a wind from a massive stel- 
lar cluster situated close to the Galactic centre. The wind 
travels through a shock and then forms a nozzle aligned 
with the rotation axis of the Galaxy, which will is to the 
gravitational potential given by Eq. (jlSp . and the environ- 
mental pressure stratification given by Eqs. (|13m4p . 

We fix the environmental pressure stratification and the 
gravitational potential as described in the previous section 
by assuming a temperature of 10^ K and a density on the 
galactic plane of 0.5 cm~'^ for the hot halo. The remain- 
ing free parameters of the model are then the cluster wind 
velocity v^, its mass-loss rate M, and metallicity Z. The 
initial flow velocity, sound speed, and density of the flow 
are functions of only the wind velocity Vw and the pressure 
of the hot environmental gas close to the galactic plane (see 
Eqs. EEl)- The metallicity Z is important for determining 
the cooling rate. The mass-loss rate is only important for 
determining the radius of the nozzle flow cross-section (see 
Eq. [5|), so that we choose a single mass-loss rate M=10~^ 
M0 yr~^ for all models. 

This is a reasonable mass-loss rate in the following sense. 
To reach a height of ^ 10 kpc above the galactic plane, the 
nozzle flow must exist for a time ~ 10 kpc/w^, '^10 Myr 
(where Vw 1000 km s~^ approximate cluster wind veloc- 
ity). This timescale coincides with the lifetime of the stellar 
population of a SSC (Leitherer & Heckman 1995). During 
this lifetime, a cluster of ~ 10000 massive stars (with a to- 
tal mass of ^ 10^ M©) will have lost ~ 10% of its initial 
mass. If we consider a much smaller cluster such as, e.g., 
the Arches cluster near the Galactic centre (which has only 
~ 100 massive stars), it would be necessary to consider a 
much lower mass-loss rate 10~^ Mq yr""'^). However, the 
mass-loss rate from the cluster enters only in the calcula- 
tion of the cylindrical radius of the nozzle flow (see Eq. 8 
and Fig. 2), and the dynamical characteristics of the flow 
do not depend on the value of the mass loss. 

We should note that a cluster of lO"* — 10^ massive stars 
will feed a ~ 10^ — 10^ Mq nozzle flow, which in princi- 
ple could produce the mass observed in an HVC complex 
(see Wakker & van Woerden 1997). A cluster with 100- 
1000 massive stars will produce a ~ lO'^ — 10^ Mq nozzle 
flow, corresponding to the typical mass of an IVC (see de 
Gouveia Dal Pino et al. 2009). 

We now carry out numerical integrations of Eqs. ([T]) and 
([3]) with the initial conditions given by Eqs. for dif- 

ferent values of v^^, and Z. We use the cooling function of 
Raymond, Cox & Smith (1976), who tabulated the cooling 
rate as a function of the temperature T and metallicity Z 
of the gas (a density dependence oc p^, appropriate to a 
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Fig. 2. The radius of the nozzle flow cross-section as a func- 
tion of height X above the galactic plane. The soHd line 
shows the Z — 1 model and the dashed line shows Z ~ 5 Zq 
model. 



low-density regime cooling, is assumed). Figures 2-5 show 
the results obtained for models with Vw — 700 km s^^ and 
Z = 1, and 5 solar metallicities. Figure 2 shows the radius of 
the nozzle flow cross-section as a function of height x above 
the galactic plane. The Z = 1 model (solid line) shows a 
radius that decreases for x < 0.10 kpc, and then increases 
monotonically at greater heights. This behaviour is similar 
to that found in adiabatic nozzles, in which the nozzle ra- 
dius decreases in the region of subsonic flow, and increases 
with distance after the sonic transition. The radius of the 
nozzle flow diverges at Xm ~ 15 kpc. This divergence was 
not found in previous nozzle solutions (see, e.g., Blandford 
& Rees 1974 and Raga & Canto 1989), because they did 
not include the effect of the gravitational force acting di- 
rectly on the material within the nozzle flow. From Figure 
2, we see that the solution for Z — 5 (with stronger radia- 
tive cooling) shows a sudden collapse of the nozzle radius at 
X « 40 kpc. This collapse is caused by the strong radiative 
cooling at temperatures below 10^ K, and is qualitatively 
similar to the results found for cooling nozzle flows by Raga 
& Canto (1989). Figure 3 shows the velocity along the noz- 
zle flow for the two models. The flow velocity flrst grows 
with increasing x, reaches a maximum, and then collapses 
to zero at the point where the nozzle radius diverges (see 
figure 2). We show in Figure 4 the density as a function of 
position along the nozzle flow. For Z = 1, a monotonically 
decresing density vs. x relationship is found. For Z = 5, 
a sharp increase in the density is found at the position in 
which the nozzle radius collapses (see Fig. 2). 

Finally, in Figure 5 we show the temperature structure 
of the nozzle flow. The Z — 5 model shows a sudden collapse 
at X « 40 kpc, which is absent in the Z — 1 model. To 
explore of parameter space, we have run models with v^, — 
500-800 km s-i (in Av^ = SOkms"^ steps) and Z = 1-10 
(in AZ = 1 steps). From these models, we computed the 
maximum height Xm reached by the nozzle flow. The values 
obtained for function of Vw and Z are shown in 

Figure 6. 

From this flgure, we see that flows with higher wind ve- 
locity reach larger heights Xm ■ Higher values of also 
obtained for decreasing metallicities. The radiative nozzle 
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flows inject cool material at Xm, which will then fall back 
onto the galactic plane in the form of dense, neutral clumps. 
When these clumps have reached a height xq, they will have 
the infall velocities given in Figure 1. 



lOL 
0.01 
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Fig. 3. The gas velocity as a function of height x above the 
galactic plane. The lines have the same description as Fig. 
2. 




Fig. 4. The numerical gas density as a function of height x 
above the galactic plane. The lines have the same descrip- 
tion of Fig. 2. 




Fig. 5. The gas temperature as a function of height x above 
the galactic plane. The lines have the same description as 
Fig. 2. 



6. Discussion: low and high velocity HI clouds 

For the adiabatic case, we obtain nozzle flow solutions that 
extend out to infinity for cluster wind velicities of > 
Vesc, where Vgsc is the escape velocity from the position of 
the stellar cluster that feeds the flow (see Section 3). For 
Vm < Vesc, the flow stalls at a finite distance Xm from the 
galaxy. As the nozzle stalls, its cylindrical radius diverges, 
producing a layer of neutrally buoyant gas at a distance 
~ Xm from the wind source. If this gas cools radiatively, 
it will become negatively buoyant and eventually fall back 
onto the galactic plane, reaching velocities similar to the 
free-fall velocity (sec Fig. 1). 

Wo then integrated the equations numerically for a ra- 
diative nozzle flow. We included the cooling function of 
Raymond, Cox & Smith (1976) tabulated as a function of 
temperature for a set of different metallicities (assuming a 
"low density regime" oc density dependence). We assume 
that the nozzle flow propagates away from the central re- 
gion of our galaxy, within an isothermal halo subject to the 
Galactic potential of Espresate et al. (2002) and Allen & 
SantiUan (1991). 

We have run a grid of radiative nozzle models for cluster 
winds with velocities 500 < Vy, < 800 km and metallic- 
ities 1 < Z < 10 Zq. From these models, wc found that the 
radiative cooling lowers the height Xm at which the noz- 
zle flow stalls. Therefore, for higher material metallicities 
within the nozzle flows (resulting in higher cooling rates), 
we obtained lower values of Xm (see Fig. 6). 

In the models with Z = I Zq , as the cluster wind ve- 



locity increases from 



500 to 750 km s^^, the value 
of Xm increases from 1 to 50 kpc. For higher values of , 
the value of Xm rapidly diverges, indicating that the nozzle 
flows do not stall (and therefore never produce negatively 
buoyant material falling back onto the galaxy). 

In the models with Z = 5Zq, the same range of Xm 
values (i.e., from 1 to 50 kpc) is obtained for cluster wind 
velocities ranging from Vyj = 620 to 770 km s~^. For Vyj 
approaching 1000 km s~^, the nozzle flows never stall, re- 
gardless of the mctallicity. 

Therefore, from our models we could obtain nozzle flows 
that stall at heights ranging from 1 to 50 kpc for cluster 
wind velocities ranging from ^ 500 to 800 km s~^ (and 
solar abundances Z = 1 — 10). 

Let us now consider low velocity clouds, which arc ob- 
served at heights xq ~ 0.3—2.5 kpc and have infall velocities 
< 100 km s~^. From Fig. 1, we see that these clouds could 
only have originated in initial heights, Xm < 4. From Fig. 6. 
we see that this height corresponds to the stalling distances 
of nozzles with cluster wind velocities < 620 km s~^ for 
Z = IZq and < 720 km s"! for Z = 10 Zq. High ve- 
locity clouds are observed at heights xq > 12 kpc and have 
infall velocities 100 200 km s~^. From Fig. 1, we see that 
clouds with infall velocities of 150 km s^^ at xo = 15 kpc 
(these being representative values for high velocity clouds) 
have originated at initial heights Xm ~ 26 kpc. From Fig. 6, 
we see that this height corresponds to the stalling heights 
of nozzles with = 720 km s"""^ (for Z = 1 Zq) up to 
Vy, = 780 km s-i (for Z = 10 Zq). 
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Fig. 6. Isocontours of maximum height Xm (in kpc) reached 
by the nozzle flow as function of the metalhcity Z and clus- 
ter wind Vyj. 

In this way, wc see that a sequence of nozzle models 
of increasing with 500 < < 800 km s~^ naturally 
explains the range of infall velocities observed from low to 
high velocity HI clouds. Nozzle flows produced by cluster 
winds with higher do not stall, and therefore will not 
produce large amounts of negatively buoyant material at 
high altitudes that could form the observed HI clouds. 

7. Conclusions 

We have presented a de Laval nozzle model for galactic 
fountain flows. In this flow, the wind emitted (in our case) 
by a stellar cluster is thcrmalized in a wind shock, and is 
then reaccelerated out of the galactic plane into a coUi- 
mated flow as a result of the pressure stratifiction of the 
surrounding environment. 

The equations for this kind of flow were flrst studied 
in the context of the production of jet flows from AGN by 
Blandford & Rees (1974), and in the context of jets from 
young stars by Konigl (1982). Raga & Canto (1989) studied 
the effect of radiative cooling within the nozzle flow (also 
in the context of jets from young stars). 

We have written the equations for a nozzle flow, includ- 
ing both a radiative cooling term and the action of the 
gravitational force on the material within the nozzle flow. 
This latter effect is caused because the gravitational force 
(in our case, because of the Galactic potential) producing 
the environmental stratification also acts directly on the 
material flowing along the nozzle. This effect has not been 
included in any of previously published astrophysical nozzle 
flow models. 

Prom numerical integrations of the radiative nozzle flow 
equations, we found that cluster winds with 600 < t;„, < 
750 km s~^ produce nozzles that stall at heights 3 < Xm < 
50 kpc above the galactic centre (sec Fig. 6). The negatively 
buoyant clumps that become separated from the stalled 
flows fall back towards the galactic plane at velocities < 
200 km s~^ (see Fig. 1). This velocity range covers the infall 
velocities found for both low and high velocity HI clouds. 

From this study, wc conclude that the high and low 
velocity clouds observed to be falling onto the Galactic disk 



could indeed be produced by winds from massive stellar 
clusters close to the Galactic centre. The condition for these 
clusters to produce the observed clouds is that the cluster 
wind velocity should be lower than ^ 750 km s^^, with 
a somewhat higher or lower velocity limit depending on the 
metallicity Z of the gas in the cluster winds (see Fig. 6). 
Clusters such as the Arches (Stevens & Hartwell 2003) and 
the Quintuplet clusters (Rockefeller et al. 2005) are clear 
candidates for driving this type of nozzle flow. 

The Arches and Quintuplet clusters eject masses of the 
order of the mass in an IVC (see Section 5). Therefore, 
several clusters such as these axe necessary to produce an 
appreciable fraction of all observed IVC's. A population 
of this clusters might be present in the Galactic centre 
(with several still undetected clusters). On the other hand, 
a complex of HVCs has a combined mass of ~ 10^ — 10^ 
Mq (sec Section 1, 5). To eject these masses from the 
Galaxy, we would need considerably more massive SSC's 
(with ^ 10* — 10'^ massive stars). It appears unlikely that if 
these clusters were present around the Galactic centre they 
would not have yet been detected. Therefore, it appears 
more likely that these clusters (if they exist) are located 
far away from the Galactic centre. It would therefore be 
interesting to explore models of nozzle flows produced by 
clusters at some distance from the Galactic centre. 

We should mention the work of Wiinsch et al. (2008). 
These authors showed that the presence of thermal insta- 
bilites inside the wind from a very massive SSC (with more 
than 8 x 10® solar masses) results in a substantial decrease 
in the total mass-loss rate of the cluster wind. This effect 
would serve to moderate the strength of fountain flows pro- 
duced by these clusters. 

To proceed with this study, it will be necessary to re- 
move at least some of the simplifying assumptions of our 
present model. For example, we have assumed that the 
nozzle flow has a uniform velocity across the nozzle cross- 
section. The removal of this assumption would probably 
lead to flows in which an outer, slower envelope could stall, 
while a higher velocity core (moving more repidly than the 
escape velocity) could continue flowing out to large dis- 
tances from the Galaxy. In this a model, cluster winds with 
higher values of might be able to produce nozzles that 
form negatively buoyant, infalling clumps. 

Inn our work, we have also assumed that the cluster pro- 
ducing the nozzle flow is located close to the Galactic cen- 
tre. It would be interesting to study off-centre cluster winds, 
which should produce nozzle flows with curved paths, be- 
cause HVCs could indeed be formed by massive clusters at 
some distance from the Galactic centre (see above). 

Finally, it would also be interesting to study the spatial 
distribution of the infalling clumps (produced by a nozzle 
flow) on the galactic plane. To address these issues, it will 
probably be necessary to carry out full, 3D gasdynamic 
simulations similar to those of de Gouveia Dal Pino et al. 
(2009, who considered the production of high and low ve- 
locity clouds by SN explosions). 
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